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Oh, Abstract 

The problem of the long-period gravity modes of the Earth outer fluid core (the 
core undertones) is investigated using either the subseismic or the anelastic approx- 
imation. These two approximations aim at filtering out the uninteresting short- 
■ period seismic or acoustic oscillations while taking into account the density vari- 

es : ations across the core. However, they differ on the form of the equation of mass 

conservation since the density perturbations do not contribute to the mass balance 
in the anelastic case. Here we show that these two approximations lead to the almost 
same results because of the weakness of the core stratification due to the convec- 
tive mixing. The anelastic approximation should be however preferred because it 
is simpler, mathematically self-consistent and is the only one which can be applied 
to problems with a stronger stratification such as the ones encountered in radiative 
q \ zones of stars. 
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1 Introduction 



To study the long-periods oscillations of the Earth outer core, geophysicists 
commonly use the so-called "subseismic approximation" . This approximation 
has been proposed by (Smylie and Rochester (1981|) to reduce the order of the 



non-rotating governing equations from fourth to second. It aims at filtering 
out the high-frequency seismic P-waves from dynamics while taking into ac- 
count the density variations across the core. Using the classical Boussinesq 
approximation for this problem is indeed quite satisfactory since, although 
the acoustic-like waves are filtered out, the equilibrium density is assumed to 
be constant ( [Spiegel and Veronis, 1960] ). 



A similar approximation, the so-called "anelastic approximation", is used in 
meteorology and astrophysics to describe compressible convection. This ap- 
proximation has first been derived by Patchelor (1953|) and |Ogura and Phillip 



(1962) to study the dry convection of the Earth troposphere. Their aim was 
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to eliminate the sound waves from the hydrodynamic equations since these 
rapid waves impose very small time steps in the numerical integration. In 
addition, they wanted to take into account the density variations across the 
atmosphere so the Boussinesq approximation was still too restrictive. Quite 
naturally, this anelastic approximation has been used in astrophysics to the 
convective zones of stars (see e.g. |Gilman and Glatzmaier, 1981|) . This approx- 
imation has also been used to describe the long-period oscillations propagating 
in the Earth atmosphere and star interiors. For example, |Durran (1989Q cal- 
culated the low-frequency oscillations of an isothermal atmosphere using the 
anelastic equations of |Ogura and Phillips (1962| ). In the same way, pintran 



and Rieutord (2000Q computed the oscillations of rapidly rotating stars. In 
order to decrease the size of the numerical rotating problem, a good solution 
indeed consists to filter out acoustic quantities while taking into account the 
large density variations encountered in radiative zones of stars. 

The subseismic approximation has also been used in the field of stellar oscil- 
lations to describe the long-period gravity modes propagating in non-rotating 
stars. |De Boeck et al. (1992| ) studied these modes by using the same subseismic 



set of equations as in |Smylie and Rochester (1981 ). More recently, pintran 



and Rieutord (2001|) proposed a comparative study of the anelastic and sub- 
seismic approximations applied to this problem (hereafter referred to as Paper 
I). There we showed that the subseismic approximation is not in fact a consis- 
tent approximation for the equations of motion. This result has been asserted 
by comparing the anelastic and subseismic eigenfrequencies with the complete 
ones in the case of the homogeneous star and the polytrope n = 3. Analytic 
solutions have been derived for the homogeneous case and the anelastic eigen- 
frequencies appear to be about twenty times more precise than the subseismic 
ones. In the same way, the computed anelastic eigenfrequencies are of about 
five times more precise than the subseismic ones in the polytropic case. We 
therefore concluded on the better efficiency of the anelastic approximation as 
far as stellar low-frequency oscillations are concerned. 

The aim of this paper is to compare now both approximations in a geophysical 
context, namely that of the long-period oscillations of the Earth outer core. 
The subseismic approximation has been already tested in core dynamics and 
the resulting errors appear to be always at most (9(1%) fllSmylie et al., 19"H3| ; 
Crossley et al., 1991| ; |Crossley and Rochester, 1992| ). No comparative study 



with the anelastic approximation has been however made for this problem. In 
particular, the real benefit of the using of the subseismic approximation instead 
of the anelastic one is not clear. In fact, we will show that both approximations 
lead to the same results for the core undertones. It means that, contrary 
to stars, the inconsistency of the subseismic approximation does not induce 
unfortunate consequences in this case. Using the anelastic approximation is 
however preferable because it is self-consistent, simpler and can be successfully 
extended to more stratified astrophysical problems. 
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The plan of the paper is the following: after setting out the basic equations 
governing the adiabatic oscillations of a non-rotating fluid (Section 2.1), we 
recall the main features of the subseismic and anelastic approximations (Sec- 
tion 2.2) and give some details on the numerical method used to solve the 
oscillation equations (Section 2.3). The results of our calculations are then 
presented with the comparisons of the subseismic and anelastic eigenperiods 
with the complete ones (Section 3). We conclude with some outlooks of our 
results, with in particular their immediate implication to the more difficult 
rotating problem (Section 4). 



2 The formalism 



2.1 The equilibrium model and the complete oscillation equations 



The identification by seismology of stable stratification zones in the Earth 
outer fluid core is still today controversial. Two stable stratified zones are in 
fact expected: (i) one zone may be located at its bottom (ICB) due to the 
release of light elements by the crystallization of the inner core ( Bouriau and 



Foupinct, 1991| ); (ii) a second zone may be located at its top (CMB) due to 



a thermal buoyancy flux coming from the mantle ([Lister and Buffet, 1998 



From observations of normal mode eigenperiods, [Masters (1979| ) constrained 
the minimum value of the core Brunt- Vaisala period to be 6h, a value con- 
sistent with the core mixing due to convection flCrossley, 1984| ; |Smylie et al. 



1984! ). Thus many authors flCrosslcy and Rochester, 1980 ; CJrossley et al., 199T| 



Urossley and Rochester, 1992| ) chose to take a uniform Brunt- Vaisala period 



of 6h in their studies of core dynamics and this model has also been used here. 

Assuming a time- dependence of the form exp(zcrt) and neglecting rotation, the 
adiabatic oscillations of a self-gravitating fluid around its reference state obey 
the following linearized equations 



\Po J \ Po9o 




r- 



— + cgdivf-0 O £r = O, (2) 
Po 

( p< N 2 \ 

A0' = AnGpo — + — £ r . (3) 
\P0c5 g J 



where the variables P' and <p' are respectively the Eulerian perturbations of 
pressure and gravitational potential and £ is the Lagrangian displacement. 
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The radial profiles po,<?o and Cq respectively denotes the density, gravity and 
square of seismic P-wave velocity and are given by the stable 6h core model 
of |Crossley and Rochester (1992| ) with 



r = -*r^ =fel = 8.461 x 10- (4) 



As in |Urossley and Rochester (1980| ), we assume rigid boundaries both at the 
ICB (r = a) and CMB (r = b) since eigenfrequencies computed using such 
conditions are in excellent agreement with those computed for the whole Earth 
(see also Crosslcy et al., 1991 ). In addition, the gravitational potential should 
be continuous across these surfaces which leads, for a mode of degree £, to the 
following boundary conditions 



£ r = — (f)' = at r = a, 

dr r 



d(j>' iU + 1) „ 
c , r = -r- + V = at r = b. 



(5) 



dr 



2.2 The inconsistency of the subseismic approximation 



The anelastic and subseismic approximations are based on the same assump- 
tion which is that the Eulerian fluctuations P' do not contribute to the La- 
grangian ones SP, that is 

dP P' 
6 P = P > + = P > _ pmir ~ _ po O £ r if ( 6 ) 

dr p g 



Quite surprisingly, this common basic assumption does not lead to the same 
approximated form of the equation of mass conservation since 



Anelastic: div^* + , ^° £ r = 0, Subseismic: div<f — ^-^ r = 0. (7) 
dr Co 



We showed in I that the subseismic form of this equation is in fact incompatible 
with the basic assumption (6) since the neglect of P' in 8P also implies the 
neglect of p' in 5p, i.e. 

D ,^^o t , dpo e x „ dP dp 

P < —ir => P < ~rCr or 5P ~ — £ r => 5p ~ — f r . 
dr dr dr dr 
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As a consequence, the only consistent form of the equation of mass conserva- 
tion is the anelastic one according to 

p' + ^ r + p div£ = 0^div£ + ^^ r = or div(p £) = 0. 



The density variations must be neglected in the mass balance if (6) is assumed. 

An other important result derived in I concerns the uncoupling of the gravity 
perturbations <p'. Assuming (6), we indeed deduce the following system from 
Eqs. (1), (7) and boundary conditions (5) 



<7 2 V x £ = V x (N 2 £ r e r ) } 
<divf+^£ r = or divf-^ P = 0, ( 8 ) 

(XT Cq 
trip) = 6(6) = 0. 

V 



This second-order set of equations valid for £ and a 2 clearly forms a well-posed 
eigenvalue problem so that the subseismic and anelastic eigenfrequencies can 
be computed without the taking into account of the 4>' perturbations. 



2. 3 Numerics 



As mentioned in Section 2.1, we choose the same equilibrium model as in 
Urossley and Rochester (1992| ), that is a stable core with a Brunt- Vaisala pe- 
riod of 6hr. This constant period gives us a natural time scale for our problem 
whereas the CMB radius b will be the length scale. Then the dimensionless 
perturbed fields are expanded onto spherical harmonics with, for example, 

+oo pi +oo 

£(r, = £ u e (r)Y e e r + v t (r)VY t , — (r, = £ Pi (r)Y £ , 

£=0 PO 1=0 



where Yt(9, ip) denotes the normalized spherical harmonic of degree I with a 
zero azimuth m (m = and no toroidal component needs to be taken into 
account since rotation is neglected; see e.g |Rieutord, 1991| ). As an example, 
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the dimensionless form of the approximated set (8) reads 



2 / dv e s 
A \u e -v e - x— 

\ (JjUj , 



x— - + 2 + Px)u t 
dx 



ui{rj) = u e (l) 



0. 



+ l)v t = 0, 



a dln P0 a 9o\ 

(3 = — — or (3 = — J , 
, dx eg I 



where rj = a/b and A 2 = a 2 /N 2 . This system may be formally written as 



M A ^£n = X 2 in M B ^i n with 



where and M. b denote two differential operators and \I> £ n is the eigenvec- 
tor associated with the eigenvalue Xj n of order n. This differential eigenvalue 
problem is discretized on the Gauss-Lobatto grid associated with Chebyshev's 
polynomials. Eigenvalues and eigenvectors are then computed with the it- 
erative Arnoldi-Chebyshev solver already used to calculate the oscillations of 
rotating stars (for more details on numerics see pintrans et al., 19991 pintran 
and Rieutord, 200q ). 



3 Results 



In order to compare both approximations, we first calculated the eigenfrequen- 
cies of the complete problem consisting of Eqs. (1-3) plus boundary conditions 
(5). Then, we compared them with their anelastic and subseismic counterparts 
computed from (8). Table 1 shows the obtained periods (in hours) for some of 
the twenty first undertones of degree I — 2, the periods being deduced from 
the dimensionless eigenvalues \\ n using I} n (hr) = 6/ \j\\ n - 

The first point is that the lower undertone period T\ ~ 9.6705 hr well agrees 
with that computed by CJrosslcy and Rochester (1992| ) who found T\ ~ 9.6741 



hr. The slight difference may be explained both by our different boundary 
conditions and the eigenvalue solver since these authors did not assume rigid 
boundaries and used a shooting method to compute eigenvalues. 

The second point is that the anelastic and subseismic approximated periods 
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Table 1 

Periods (in hours) of some of the twenty first I = 2 undertones of the stable 6hr core, 
with their corresponding anelastic and subseismic counterparts and the involved 
errors in percent. 



n 


Complete 


Anelastic 


Errors (%) 


Subseismic 


Errors (%) 


1 


9.6705371 


9.6302857 


0.4162 


9.6328597 


0.3896 


3 


23.031630 


23.012466 


0.0832 


23.013470 


0.0788 


5 


37.465939 


37.453979 


0.0319 


37.454590 


0.0303 


10 


74.139875 


74.134154 


0.0077 


74.134462 


0.0073 


15 


110.98888 


110.98481 


0.0036 


110.98502 


0.0035 


20 


147.88150 


147.87863 


0.0002 


147.87878 


0.0002 



are almost the same. This can be easily understood from the different forms 
(7) of the equation of mass conservation where the £ r -term is either equal 
to (dhipo/dr)£ r or — (<7o/ c o)£r- Because the Brunt- Vaisala frequency is weak, 
these two terms are almost equal; indeed from (4), we have 



go d\np 







dr 



AT* 

go 



8.5 x 10" 8 s- 2 
780 cm/s 2 



-1(T 5 km 



-i 



The difference in accuracy between both approximations thus depends on the 
stratification strength. As the Earth outer core is almost adiabatically strat- 
ified due to convection, the subseismic and anelastic approximations give the 
same results, which are very close to the exact ones as shown by Table 1. As 
expected, the relative errors decrease with increasing order of the mode. The 
accuracy level of both approximations is remarkable since, except for the first 
undertone, relative errors are less than 0.1% ! In fact, this very good agreement 
is still related to the weak stratification. When applying the anelastic approxi- 
mation, one should indeed be careful about the required physical assumptions. 
Batchclor (1953 ) and Pgura and Phillips (1962p derived their original anelastic 



set of equations using two assumptions: 

(1) the reference state around which oscillations occur is almost adiabatically 
stratified. 

(2) the time scale of any disturbance is similar to that of gravity waves. 

Among these two basic assumptions, the second one is the easiest to fulfil. 
In this paper, we chose the constant Brunt- Vaisala period of our equilibrium 
model as the time scale. It automatically satisfies the second assumption since 
the computation is restricted to disturbances such as d/dt ~ N. 

The first assumption is more severe. Assuming an almost adiabatic state means 
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that the equilibrium potential temperature 6q is also nearly constant; i.e. 
s = 56 q /© <C 1 where 59 is the jump of the equilibrium potential tempera- 
ture accross the layer and G is the reference mean value. Pgura and Phillip^ 
(1962|) expanded all dependent variables as a power serie of e and obtained the 
anelastic equations by collecting the first-order terms. The neglected terms are 
then 0{e 2 ) and give the accuracy of the development. This small dimensionless 
parameter e may be related to the Brunt- Vaisala frequency by 

^ 2 _ din go ^ go 56q ^ 

dz d d ' 



where d denotes the layer thickness. In our case, we have 

d „ r9 2000 km _ _ n „ 9 nM 

e ~ — N 2 ~ 9 8.5 x 10~ 8 s~ 2 ~ 2%. 

£o 780 cm/s 



Hence the mean error made when using the anelastic approximation on the 
core undertones is about 0(e 2 ) ~ 0.04% which agrees well with the results 
from Table 1. The long-period Earth core oscillations can therefore be seen as 
a good case for the anelastic approximation and, despite it inconsistency, also 
the subseismic approximation since the almost neutral stratification prevents 
unfortunate consequences in this case. On the contrary, as the stratification 
increases like for instance in stellar radiative zones, these two approximations 
disagree more and more and the anelastic results become more accurate as 
shown in I. 



4 Conclusion 



We have computed the complete, anelastic and subseismic approximated eigen- 
periods of a stable 6h Earth's outer core. After setting out the complete equa- 
tions governing the adiabatic oscillations of a non-rotating and self-gravitating 
fluid, we presented the basic properties of the anelastic and subseismic approx- 
imations. 

Both approximations assume that the Eulerian pressure perturbations P' do 
not contribute to the Lagrangian ones SP, that is only the fluctuations stem- 
ming from the equilibrium pressure gradient are retained. We then recalled 
that, unlike the anelastic case where the perturbed density field has no influ- 
ence on the mass balance, the subseismic equation of mass conservation is not 
compatible with this assumption (this result being demonstrated in I). As a 
consequence, the subseismic approximation is not a consistent approximation 
of the equations of motion. 
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We also recalled that the gravity perturbations decouple from the motion in 
both cases. It refutes the common belief in the respect of self-gravity by the 
subseismic set of equations; i.e. the only difference with the anelastic approxi- 
mation lies in the equation of mass conservation. Then, two main results have 
been enlightened by our computations: 

• for the Earth's fluid outer core, anelastic and subseismic eigenperiods are 
very close to each other. We explained this result by the very weak strat- 
ification of the fluid. The differences between the anelastic and subseismic 
periods indeed come from the stratification strength. Because of convection, 
the Earth core is almost adiabatically stratified and both approximations 
give the same results. 

• the anelastic relative errors are very small (<,0.1%) for core undertones. It 
is not obvious at first sight since the anelastic set of equations has been 
initially derived for the atmospheric dry convection. We explained this very 
good agreement by recalling the required physical assumptions implied by 
the use of the anelastic approximation: (i) the reference state should be al- 
most adiabatically stratified; (ii) the time scale of any disturbance should be 
similar to that of gravity waves. The core undertones clearly fulfil very well 
these two basic assumptions and the anelastic eigenperiods are thus accu- 
rate. Because of the weakness of the Brunt- Vaisala frequency, the subseismic 
results are the same as the anelastic ones and both approximations can be 
in fact applied to this problem. However, it is more natural to employ the 
anelastic approximation since it is simpler, mathematically self-consistent 
and it is the only one which can be successfully applied to problems with a 
stronger stratification such as the long-period oscillations of stars. 

The above results have straightforward consequences as far as geophysical ap- 
plications are concerned. The main application concerns the rotating problem, 
that is to calculate the eigenperiods of the rotating outer fluid core, pintrans 
et al. (1999| ) already calculated the gravito-inertial waves propagating in a 
rapidly rotating stratified shell which had the same aspect ratio than the 
Earth outer core (i.e. rj = a/b ~ 0.35). To decrease the size of the numeri- 
cal problem, this work has been done under the Boussinesq approximation so 
that the density variations across the outer core have been neglected. Hence 
an improvement of this problem is now possible by taking the more complete 
equilibrium model used in this paper by the means of the anelastic approxi- 
mation. 
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